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Abstract: QCD was shown to have a nontrivial vacuum structure due to the topol- ^ ^ 

ogy of the 6 = 9 + 27m parameter. As a result of this nontrivial topology, in the 1^^^ 
large limit, quasi-stable QCD domain walls appear, characterized by a transition CD 
in the singlet rj' field. We discuss the physics of these QCD domain walls as well CD 

as related axion domain walls and we present a new type of axion wall which also | \ 

contains an rj' transition. We argue that these domain walls are topologically stable >. ^ 

in the limit Nc oo and classically stable for large but finite Nc, however, they can 
decay through a tunneling process. We argue that the qualitative features of these 
QCD domain walls — namely their classical stability — persist to the realistic case ' ' 

of A'c = 3 and that it is at least possible that their lifetime could be macroscopically 
large. If it is, then QCD domain walls could play an important role in the evolution 
of early universe and may be detectable in energetic collisions such as those at the 
Relativistic Heavy Ion Collider (RHIC). 
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1. Introduction 

Colour confinement, spontaneous breaking of chiral symmetry, the U(l) problem, 6 
dependence, and the classification of vacuum states are some of the most interesting 
topics in QCD. Unfortunately, the progress in our understanding of them is extremely 
slow. At the end of the 1970s A. M. Polyakov demonstrated colour confinement in 
3-dimensional QED (QED3): this was the first example in which nontrivial dynamics 
of the ground state played a key role. Many papers were written regarding the ground 
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state structure of gauge theories in the strong couphng regime, but there were many 
unanswered questions. Almost 20 years passed before the next important piece of the 
puzzle was solved p|, |^ . Seiberg and Witten demonstrated that confinement occurs 
in supersymmetric (SUSY) QCD4 due to the condensation of monopoles: a similar 
mechanism was suggested many years ago by 't Hooft and Mandelstam (see IQ for a 
review). Furthermore, condensation of dyons together with oblique confinement for 
nonzero vacuum angle 9 was also discovered in SUSY models (this phenomenon 
was also argued to take place in ordinary QCD. See [Q). 

In addition to providing solid demonstration of earlier ideas, the recent progress 
in SUSY models has introduced many new phenomena, such as the existence of rich 
vacuum state structure and the existence of domain walls 0: topologically stable 
interpolations connecting the same vacuum state. The same conclusion was reached 
by Witten in |^ based on a D-brane construction in the limit of large N^.. In fact, 
one can see that in both approaches, the number of states in the vacuum family is N^- 
Motivated by this development in SUSY gauge theories D-brane construction, we ask 
how this applies to QCD. The issue of classifying the vacuum states for a given param- 
eter 9 as well as the phenomenological consequences (domain walls and their decay 
etc. ) is the subject of this work. In separate publications 0] we shall apply these re- 
sults to analyze the physics after the QCD phase transition in evolution of early Uni- 
verse. In particular, we discuss the possibility of generation of primordial magnetic 
fields due to the existence of short-lived QCD domain walls discussed in this paper. 

The starting point of our analysis is an effective Lagrangian approach. Expe- 
rience with SUSY models demonstrates that the effective Lagrangian approach is a 
very effective tool for the analysis of large distance dynamics in the strong coupling 
regime. There are two different definitions of an effective Lagrangian in quantum 
field theory. One of them is the Wilsonian effective Lagrangian describing the low 
energy dynamics of the lightest particles in the theory. In QCD, this is implemented 
by effective chiral Lagrangians for the pseudoscalar mesons. Another type of effective 
Lagrangian is defined by taking the Legendre transform of the generating functional 
for connected Green's functions to obtain an effective potential. This object is useful 
in addressing questions about the vacuum structure of the theory in terms of vac- 
uum expectation values (VEVs) of composite operators — these VEVs minimize the 
effective action. The latter approach is well suited for studying the dependence of 
the vacuum state on external parameters such as the light quark masses or the vac- 
uum angle 9. However, it is not, for example, useful for studying S'-matrix elements 
because the kinetic term cannot be recovered through this approach. The utility of 
the second approach for gauge theories had been recognized long ago for supersym- 
metric models, where the anomalous effective potential was found for both the pure 
supersymmetric gluodjTiamics [|I^ and supersymmetric QCD (SQCD) [|1T| models. 
Properties of the vacuum structure in the SUSY models were correctly understood 
only after analyzing this kind of effective potential. 



This paper contains many of the details mentioned in the letter and is orga- 
nized as follows: 



Section Here we review the properties of the QCD effective Lagrangian [|T2|, 
13|1 which is a generalization of the Di Vecchia-Veneziano-Witten (WW) La- 



grangian [jT4|, |l^] to include terms sub leading in l/Nf. as well as to account for 
a constraint due to the quantization topological charge^. One should empha- 
size from the very beginning that the specific form of the effective potential 
used in this paper is not critical for the present analysis: only the topological 
structure and winding 0^6 + 27in — a consequence of the topological charge 
quantization — is essential. 

Section ^: Here we review the idea of topological charge conservation and argue 
that, in the limit Nc — > oo, the QCD effective Lagrangian admits stable domain 
walls. We then discuss how these walls are metastable on the quantum level 
due to a tunnelling decay mode allowed when we consider heavier degrees of 
freedom to the next order in l/Nc- 

Section ^: Here we present the classical domain wall solutions for QCD domain 
walls both with and without a dynamical axion field. 

Section ^: Here we show that the QCD domain walls are not stable on the quantum 
level due to a tunnelling phenomena as discussed in Section ^. In this section 
we also discuss the validity of the large Nc approximation when Nc = 3 and 
argue that the realistic case is qualitatively the same as the large Nc limit. We 
estimate the lifetime of the QCD domain walls when A^^^ = 3 and show that 
it is potentially macroscopically large. Thus, these domain walls may play an 
important role in various physical processes. 

Section Bl: This is our conclusion. 



2. Effective Lagrangian and 9 dependence in QCD 



Our analysis begins with the effective low energy QCD action derived in |]T2|, |T^], 
which allows the 6'-dependence of the ground state to be analyzed. Within this 
approach, the pseudo-Goldstone fields and rj' field are described by the unitary ma- 
trix Uij, which correspond to the 75 phases of the chiral condensate: (\E'^\Ef'jj) = 
-\(^l'^r)\U,j with 



U = exp 



iV2— hi- 



V' 



uu^ = 1 



(2.1) 



^Such a generalization was also motivated by SUSY consideration ITfl], see also H for a review. 



where A*^ are the Gell-Mann matrices of SU(iVj-), tt" is the pseudoscalar octet, and 
/tt = 133 MeV. In terms of U the low-energy effective potential is given by [ll2| , 0: 



^"'^ f^rT./,^rr XT N T^r. / 27r/ + i log Dct t/ 




1 v-^ V 

lim — log V exp< — Tr(Mf/ + H.c.) + VEcos , 



All dimensional parameters in this potential are expressed in terms of the QCD 
vacuum condensates, and are well known: M = diag(m^|(^*\l/*)|); the constant 
E is related to the QCD gluon condensate E = i^G"^), where numerically 3b = 
llNc — 2Nf; the quark condensate (^*\E'*) ~ —(240 MeV)^, and the gluon condensate 
{^G^) - 1.2 X 10-^GeV\ 

It is possible to argue that Equation ( p.2| ) represents the anomalous effective 
Lagrangian realizing broken conformal and chiral symmetries of QCD. The arguments 
are that Equation (|2.2|): 



1. correctly reproduces the WW effective chiral Lagrangian, [^4|, [1^ in the large 
Nc limit, 

[For small values of {9 — i log Det U) , the term with / = dominates the infinite 
volume limit. Expanding the cosine (this corresponds to the expansion in 1/Nc), we 



recover exactly the WW effective potential |14, |15[] together with the constant term 



-E = —{bas/{32iT)G'ij^) required by the conformal anomaly: 



Wvvwi0,U,U^) = -E - ^TriMU + R.C.) + ^{i^'^)YMinogBetU - ef . . . , (2.3) 

where we used the fact that at large E/N^ = —{v'^)ym is the topological sus- 
ceptibility in pure YM theory. Corrections in 1/Nc stemming from Equation (|2.2|) 



constitute a new result of |12, IT 



reproduces the anomalous conformal and chiral Ward identities of QCD, 

[Let us check that the anomalous Ward Identities (WI's) in QCD are reproduced 
from Equation (p.2|). The anomalous chiral WI's are automatically satisfied with the 



substitution 9 ^ (9 — i log Det U) for any Nc, in accord with [14, [l^. Furthermore 



it can be seen that the anomalous conformal WI's of |17| for zero momentum correla- 



tion functions of the operator in the chiral limit rriq are also satisfied when 
E is chosen as above. As another important example of WI's, the topological sus- 
ceptibility in QCD near the chiral limit will be calculated from Equation (p.2|). For 
simplicity, the limit of S\J{Nf) isospin symmetry with Nj light quarks, niq <C Aqcd 



will be considered. For the vacuum energy for small 6 one obtains [T^, 13 1 



E,ac(^) = -E + mq{m)Nf cos (^^^ + 0{ml) . (2.4) 



Differentiating tliis expression twice with respect to reproduces the chiral Ward 
identities H, |l|]: 

hm i I dxe^^-(0|r{^GG(x) ^GG(O)} |0) = -^^^ (2.5) 



Other known anomalous WI's of QCD can be reproduced from Equation ( |2.2| ) in a 
similar fashion. Consequently, Equation (^]^) reproduces the anomalous conformal 
and chiral Ward identities of QCD, and in this sense passes the test for it to be the 
effective anomalous potential for QCD.] 



reproduces the known 9 dependence |14, 15]. 



[As mentioned earlier, our results are similar to those found in [14, 15|. A new 
element which was not discussed in 80's, is the procedure of summation over / in (p.2|). 
As we shall discuss in a moment, this leads to the cusp structure of the effective 
potential which seems to be an unavoidable consequence of the topological charge 
quantization.^ These singularities are analogous to the ones arising in SUSY models 
and show the non-analyticity of phases at certain values of 9. The origin of this non- 
analyticity is clear, it appears when the topological charge quantization is imposed 
explicitly at the effective Lagrangian level.] 

In general, the 9 dependence appears in the combination 9/Nc, (see Equa- 
tion ( p.2|) ) which naively does not provide the desired 27r periodicity in the physical 
observables. Equation p.2| ), however, explicitly demonstrates the 2n periodicity of 
the partition function. This seeming contradiction is resolved by noting that in the 
thermodynamic limit, V oo, only the term of lowest energy in the summation 
over / is retained for a particular value of 9. The result is that the local geometry of 
any particular 9 state, gives the illusion of 9/Nc periodicity in the observables, but 
when one considers the full topology of all the 9 states and properly switches to the 
lowest energy branch, one regains the true 9 periodicity. Of course, the values 9 and 
9 + 27rn are physically equivalent for the entire set of states, but relative transitions 
— switching branches — between different 9 states have physical significance. It is 
exactly these transitions — resulting from the non-local effects of the topology of the 
fields — that are responsible for the domain walls we discuss in this paper. 

When considering only the lightest degrees of freedom, as we do in the thermody- 
namic limit — oo of ( |2.2| ) , the effective potential acquires a cusp singularity where 
one switches from one branch of the potential to another. These cusps represent 
physical transitions in heavier degrees of freedom which have been integrated out 



^This element was not explicitly imposed in the approach of the procedure was suggested 

much later to cure some problems in SUSY models (see |16| and references therein). Analogous 
constructions were discussed for gluodynamics and QCD in O. 



to get the low-energy effective potential (|2.2|). The physical effects of these heavier 
degrees of freedom will be discussed in Sections and |5]1. The reader is referred 



to the original papers [|T2|, |T3[ for more detailed discussions of the properties of the 
effective potential (p.2|). 

Our final remark regarding ( |2.2| ). The appearance of the cosine interaction, 

's frozen): the 



cos{9/Nc), implies^ the following scenario in pure gluodynamics 
{2ky^ derivative of the vacuum energy with respect to 6, a.s 6 ■ 
solely in terms of one parameter, 1/Nc, for arbitrary k: 



0, is expressed 



QQ2k 



2k 



lld'x^{Q{x,)...Q{x2k))-^ 



2k 



[2.6) 



where, Q = ^GG. This property was seen as a consequence of Veneziano's solution 
of the U(l) problem The reason that only one factor appears in Veneziano's 
calculation is that the corresponding correlation function, ~ J Ylitid'^^iiQi^i) ' ' ' 
Q{x2k)), becomes saturated at large distances by the Veneziano ghosts whose contri- 
butions factorize exactly, and was subsequently interpreted as a manifestation of the 
6/Nc dependence in gluodynamics at small 6. However, at that time it was incorrectly 
assumed that such a dependence indicates that the periodicity in 9 is proportional 
to 2ttNc- We now know that the standard 27r periodicity in gluodynamics is restored 
by the summation over / in ( |2.2| ) such that one jumps from one branch to another 
at 6^ = TT. 

2.1 Vacuum states 

In the next section we shall discuss different types of domain walls which interpolate 
between various vacuum states, but first we should study the classification of vacuum 
states themselves. In order to do so, it is convenient to parameterize the fields U as 



U 












i4>N 



\ 



such that the potential (|2 



V 

takes the form 
-Ecos (—6 - — I 



[2.7) 



Mi cos 0j 



The minimum of this potential is determined by the following equation: 



sm 



-9--E 



E 



sm ( 



AT 



:2.8) 



[2.9) 



^As we noticed in the introduction, the specific form of the potential is not very essential for 
what follows. However, this form is very appealing for the present study because with it, we can 



describe some of the domain walls in the analytical form (see below, for example ( 4.12| ) 



At lowest order in l/N^ this equation coincides with that of |T^, |T^. For general 
values of Mi/E, it is not possible to solve Equation ( p.9|) analytically, however, in 
the realistic case e^t, <^ 1 , ~ 1 where Si = NcMi/E, the approximate solution 
can be found: 

mrfsin^ , ^ 

[m^ + + Zmumd cos u\ ^/ ^ 

m„sin6' , r>/ a 

sm 0d = — ^— — - + 0(e„, Ed) , 

[m„ + + Imurrid cos U\ ^/^ 

sin0, = 0(£„,£d). (2.10) 

This solution coincides with the one of [0, 15| to leading order in Eu,Ed- In what 
follows for the numerical estimates and for simplicity we shall use the SU(2) limit 
= ^ TTT's where the solution ( p.lO|) can be approximated as: 



(pu^O/2, (^,-0, 0<^^<7r, (2.11) 

0„ ~ (^ + 27r)/2, (pd^ {0 -27t)/2, 0, ~ 0, vr < ^ < 27r, etc. 

Once solution ( |2.11| ) is known, one can calculate the vacuum energy and topologi- 
cal charge density Q = (0|^GG|0) as a function of 9. In the limit = = m, 
{dd) = {uu) = (f'^) one has: 

Kac(^) ^ Vg=o + 2m\{^^ (l - ICOS^ 



2' 

{e\I^GG\e) = -^^^ = -m|(0|M|0)| sin ^ . (2.12) 



As expected, the 6 dependence appears only in combination with m and goes away 
in the chiral limit. One can also calculate the chiral condensate (^'};,^^) in the 6 
vacua using solution ( p.ll| ) for vacuum phases: 

(^|M|^)=cos^(0|M|0)e=o, 

{e\m^,^\9) = - sin ^ {0\m\0)e=o • (2.13) 

A remark is in order. As is well known, in thermal equilibrium and in the limit 
of infinite volume, the \6) vacuum state is a stable state for all values of 6. Thus, it 
is possible to conceive of a world with ground state \6) where 6^0. The physics of 
this world would be quite different from that of our own: In particular, P and CP 
symmetries would be strongly violated due to the non-zero value of the P and CP 
violating condensates (|2.12| ), (|2.13|) . Despite the fact that the state 16*) has a higher 



energy than |0) ( p.l2| ), it is stable because of a superselection rule: There exists no 
gauge invariant observable A in QCD that can communicate between different \6) 
states, i.e. {6'\A\6) ~ 6{6 — 6') for all gauge invariant observables A. Therefore, there 
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are no possible transitions between these states ||2T], ^ ^ and any such state \6) 
may happen to be the ground state for our world. 

On account of this superselection rule, one might ask why = is so finely tuned 
in our universe. Indeed, within standard QCD, there is no reason to prefer any partic- 
ular value of 9. This is known as the strong CP problem. One of the best solutions 
to this problem has been known for many years: introduce a spontaneously bro- 
ken symmetry (Peccei-Quinn symmetry [^). The corresponding pseudo-Goldstone 



particle — the axion ||T^, [^-[^ — behaves exactly like the parameter 6 but is 
now a dynamical field, thus we can absorb the parameter 9 by redefining the axion 
field a{x) and set ^ = 0. The axion is now dynamical and so the corresponding 
states 1^) ~ \a{x)) are no longer stable: the axion field relaxes to the true minimum 
9 ~ a{x) = (|2.12|) . The axion is included in the potential (p.2|) through the Yukawa 



interaction 

Wtt{MU + H.c.) ^ VTr(Mf/e^" + H.c.) (2.14) 

and kinetic term /^((9^a)^/4. Examining the potential ( p.2|) we see that the param- 
eter 9 can be absorbed into the U(l) phase of U which in turn can be removed 
by a redefinition of the axion field a a + 9/Nf. Although axions have not been 
detected and experiments have ruled out the possibility of the original electroweak 



scale axion p5| , P6| , there is still an allowed window with very small coupling con- 
stant fjr/fa ^ 1 emphasizing that the axion arises from a very different scale than 
the electroweak or QCD. Axions with this scale, however, are strong dark matter 



candidates (see for example The axion thus provides a way for the vac- 

uum state to relax to the lowest energy state |^ = 0) ( p.l2|) . In the following we shall 
consider two types of domain walls: Axion domain walls where 9 is the dynamical 
axion field a and QCD domain walls where 9 is the fundamental parameter of the 
theory. In the latter case, we assume that the strong CP problem has been solved 
by some means and take ^ = to be fixed. 

The effective potential ( |2.2| ) can be used to study the vacuum ground state ( |2.12| ), 



( |2.13| ) as well as the pseudo-Goldstone bosons as its lowest energy excitations. In 
particular, one could study the spectrum as well as mixing angles of the pseudo- 
Goldstone bosons by analyzing the quadratic fluctuations in the background field 
(|l2D, (|1|). We refer to the original papers [0, |3J] on the subject for details. 



but here we want to quote the following mass relationships for the rj' meson to be 
used in the following discussions: 

/X, = -j^E + —J2^,\m^^'\0)\+O{ml), (2.15) 



/ 



where i runs over the three flavours {u, d, s}. This relation is in a good agreement 
(on the level of 20%) with phenomenology. In the chiral limit this formula takes 
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especially simple form 



9 4:Nf , , 

< = 7^^' (2-16) 

J IT C 

which demonstrates that, in the chiral limit, the rj' mass is proportional to the 
gluon condensate, and is therefore related to the conformal anomaly. What is more 
important for us in this paper is that the combination on the right hand side of 
this equation exactly coincides with a combination describing the width of the QCD 
domain wall (see Equation ( [4.13D ). For this reason, the properties of the QCD domain 



walls are dominated by the rj' field. 
2.2 Large A^^ limit 



At this point, we pause to consider the large limit first discussed by 't Hooft 
In order to define the limit, we must hold g'^N^ ~ 1. From this it follows that 
tts ~ S'^ ~ Similarly, by analyzing the structure of the appropriate correlation 

functions, we find the following leading Nc dependences: 



as ~ AT-i , 6 ~ AT, , (GG) ~ A^^ , 



(2.17) 



Applying these to (p.l6|) we reproduce Witten's famous result |36 



lim m^/ ~ N-^/^ . (2.18) 



Here we have first taken the chiral limit rriq — >■ and then the large Nc limit. In 
general, when rriq ^ 0, there are 0(1) corrections arising from the second term 
in ( p. 151) , but in reality, the chiral limit is good due to the small quark masses and 



we wish to preserve the qualitative behaviour of this limit as we take Nc to be large. 
Thus, in the following discussion, we always take the chiral limit first so that, to 
leading order in rriq and Nc, the first term of ( |2.15|) dominates and ( |2.18|) holds. 



Most importantly, we see that, in the large Nc limit, the r]' becomes light, and 
thus we are justified in including its dynamics in the low-energy effective theory (|2.2| ). 
We shall return to this issue later in Section f).5\ where we will use the large Nc limit 
to regain theoretical control. 



3. Topological stability and instabilities 



The domain walls that we will discuss are examples of topological defects: classical 
solutions to the equations of motion which are stable due to the topological configu- 
rations of the fields. Examples of topological configurations abound in the literature, 
for example, Instantons, Skyrmions, Strings, Domain walls etc. [p7|, BH 
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The basic idea is that the theory contains some conserved charge density that 
is a total derivative = d^f{x). In this case, the total charge Q = Jqo d^x is 
quantized and represents some topological property of the fields such as winding or 
linking. Magnetic charge in the Georgi-Glashow model is an example (see |]3^ for 
details). 

The essential point is that the charges Q are exactly conserved quantities: the 
subspace of configurations with Q = Qiis orthogonal to the subspace with Q = Q2 7^ 
Qi. In particular, Q2 — Qi = n so there is no continuous way to vary a configuration 
from one subspace to the other and thus there is no overlap. Thus, objects with 
non-zero charges are absolutely stable: even if have an energy higher than the true 
vacuum state where Q = 0. 

Referring back to our effective theory, we shall show that there is a conserved 
topological charge associated with domain wall configurations. Physically, the states 
6 = and 6 = 27rn are identical: they represent the same state. However, as one can 
see from Equation ( |2.11| ), the solutions for the ground states corresponding to 6 = 
and 9 = 2tt are not described in the same way: ^ = corresponds to 0^ = 0^ = 
while 6 = 271 corresponds to = 27r, (p^ = 0. It clear that the physics in both these 
states is exactly the same: If we lived in one of these state and ignored the others, 
then we could assign an arbitrary phases ~ 27m for each 0^ or 0^ separately and 
independently. However, if we want to interpolate between these states to get feeling 
about both of them, the difference in phases between these states can no longer be a 
matter of choice, but rather is specified by Equation ( |2.11| ). This classification arises 
because the singlet combination (ps = really lives on a U(l) manifold which 

has the topology of a circle. The integer n is only important if we are discussing 
transitions around the U(l) circle: in this case, it is important to keep track of how 
many times the field winds around the centre. Thus, n is a topological winding 
number which plays an important role when the physics can interpolate around 
the entire U(l) manifold. We illustrate this idea in figure ^ Here, we show three 
topologically distinct paths in a two-dimensional space with an impenetrable barrier 
in the centre. The paths that wind around the barrier cannot be deformed into the 
other paths. Each path is characterized by a winding number n. 

An example of the situation described above is the well-known 1 + 1 dimensional 
sine-Gordon model defined by the Lagrangian £sg = ^{d^cpY — Acos(0). Here, the 
topological current = ^e^^^d'^cp and is related to the + 27rn symmetry of the 
ground state. The well known soliton and antisoliton (kink) solutions are absolutely 
stable objects with topological charges Q = ±1. These cannot be recovered by the 
standard methods of quantum field theory when one starts from the vacuum state 
(0) = and ignores the topology. 

In our case an analogous ground state symmetry is realized by the procedure of 
summation over / in Equation ( p.2| ), which makes symmetry 6^6 + 27m explicit. 
Thus, in 1 + 1 spatial dimensions we have the analogous stable objects. The fact 
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Figure 1: Three examples of transitions that are topologically different. Paths A and B 
wind around the hold once (n = 1) and twice (n = 2) respectively whereas path C does 
not wind (n = 0). Path C can be shrunk to a point whereas the others cannot. Each path 
is said to belong to a different homotopy class. Only paths A and B are affected by the 
topology of the space. Path C might imagine that it is living in a space with no hole. 

that we actually consider 3 + 1 dimensions means that the objects are not point- 
like solitons as they are in the sine-Gordon model, but rather, are two-dimensional 
domain walls with finite surface tension. 

3.1 Heavy degrees of freedom 

What we have said up to this point is well known. In the sine-Gordon model, 
the solitons are absolutely stable objects as can be seen by the fact that they are 
associated with the conserved current = e^,y9^0 (here, the indices run over the 
1 + 1 dimensions z and t). The conservation is trivial: = -^e^ud^d^fj) = 0. The 
corresponding topological charge is Q = / qtdz = ^ J —dz4>{z)dz = ^{4>z=~oo — 
(j)^^^^) = n_ — which is described by the winding number n = — n_. Here, 
the field is in a vacuum state (j)z=±oo = 2TTn± at infinity. Thus, we see that the charge 
is absolutely conserved and is integral. 

In our effective theory ( |2.2| ), we consider an analogous conserved current 

27rg^ = -ze^^Tr U^d^U = e^^d" ^ 0^ = e^,9>5 , (3.1) 
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where we have introduced the notation (ps = Xli 0« ~ v' that we shall use later for the 
isotopical singlet (77') field. By the same argument, we see that the two dimensional 
domain walls of the theory ( |2.2|) are absolutely stable.^ 

In the sine-Gordon model, this is the end of the story: the solitons are absolutely 



stable. In our effective Lagrangian (|2.2|), however, we have neglected the gluon 



degrees of freedom. In reality, however, the gluon degrees of freedom are not very 
heavy. Thus, we must consider these extra degrees and look at how they affect the 
charge conservation. What we find is that, when we account for the extra gluon 
degrees of freedom, the topology of the fields is no longer restricted to the U(l) 
manifold. These extra degrees of freedom allow the domain walls to continuously 
deform and to decay so that the ground state exists everywhere. 

To see how an extra degree of freedom can change the topology, consider figure ^ 
Here we have added a third dimension to show that the barrier was actually a peg of 
finite height. Now that we can move in the extra dimension, we can use this degree 
of freedom to "lift" the paths over the peg: thus, they are no longer topologically 
stable. In the QCD analogue, paths A and B represents domain walls (path C would 
is a trivial closed loop which would relax to a point representing the same vacuum 
state everywhere with no domain wall.). There is an energy cost to "lift" the path 
over the barrier and at low temperatures T <C Aqcd there is not enough energy to 
do this, so classically, the domain walls are stable. It is still possible, however, for 
the walls to overcome the barrier by tunnelling through the barrier. The tunnelling 
probability, however, could be low due to the height of the obstacle and hence the 
lifetime of the walls could be much larger than the Aq^j^ scale which one might 
naively expected for standard QCD fiuctuations. See Section |^ for details of the 
dynamics of the gluon fields. 

We should also remark that with each winding, the domain walls become more 
energetic. Walls with a large number of windings either have enough energy to 
rapidly unwind, or else separate spatially forming several domain walls of winding 
number ~ ±1. For this reason, we shall discuss in this paper only the simplest walls 
which wind once. 

To tie this picture together, consider the formerly conserved current (|3.1|) and 
effective Lagrangian (p.2|) and ask: from where do the phase fields 0j come? These 
phases arise from some sort of complex field $ = pe^'^^ . In general, we must con- 
sider the dynamics, not only of the phase 05, but of the component p = \^\. We 
assume that this field lives in some sort of Mexican-hat potential with approxi- 



■*Onc might think that, since the domain walls directly involve the 77' field, that the stability of 
the 77' particle might affect the stability of the domain walls. This is not so. Even in the effective 



theory (2^), the 77' particle can be considered as unstable decaying 77' 2j for instance. This 



instability is related to the fact that the 77' number charge is not conserved. Irrespective of this 



non-conservation, the current (3.1) is still perfectly conserved. The only way for these domain walls 



to decay is by violating the conservation of this current. This is what we consider next. 
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Figure 2: Here we show the same picture as in figure |^ except that we show the third 
dimension. Here we can see that all the paths are now homotopically equivalent. We can 
deform the paths by "lifting" them over the obstacle so that we can unwind them. If the 
paths were strings with some weight, then it would require some energy to "lift" the strings 
over the obstacle. If this energy was not available, then we would say that, classically, the 
configurations that wind around the peg are stable. Quantum mechanically, however, the 
strings could still tunnel through the peg, and so the configurations are unstable quantum 
mechanically. The probability that one string could tunnel into another configuration 
would depend on the height of the peg. 



mate symmetry 05 ^ 05 + a and minimum valley where (p) = 1. This symmetry 
is thus spontaneously broken by the vacuum expectation values of the condensate 
($) = exp(2(05)). To recover the effective Lagrangian: we "integrate" over the 
"heavy" p degree of freedom by setting p = (p) = 1 equal to its classical expectation 
value and ignoring the quantum fluctuations about this minimum. This must repro- 
duce the effective potential ( |2.2| ) with the pseudo-Goldstone field 05 (77' and vTj in 
the real theory). The appropriate Mexican-hat potential reproducing ( |2.2| ) is given 
in 

One can now consider an appropriate generalization of the current (|3.1| ) 



27rg^ = -2e^.Tr ^^8"^ = e^, (-«p9> + p'd'ct^s) (3-2) 
which reduces to ( |3.1| ) when we set p = (p) = 1. Now we have 

vrS^g^ = p e^, d'' p d''<j)s (3.3) 



which is no longer conserved as expected. This was only conserved in the effective 
theory (|2]^) because we integrated out the heavy degrees of freedom by setting p = 1 
so that d^p = 0. Thus, the decay of the domain walls is directly related to the 
dynamics of the heavy degrees of freedom. We consider this effect in Section ^ 
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The physical object responsible for this behaviour in ( |2.2|) is the gluon condensate 
E ~ (G^), which is the most essential contribution to the mass m^/ of the t]' particle.^ 



4. Domain walls 



In the rest of this paper we limit ourselves with the simplest case Nf = 2 and neglect 
the difference between /^^ and fr^' which numerically are very close to each other. To 
describe the basic structure of the QCD domain walls as well as that of axion do- 
main walls we replace the parameter 6 in Equation (|2.2| ) by a dynamical axion field 
6 — ^ Nfa = 2a (this corresponds to the so-called N = 2 axion model). We also intro- 
duce here the following dimensionless phases, 0^ describing the isotopical "singlet" 
field, and 0t describing the isotopical "triplet" field. These fields correspond to the 
dynamical rj' (singlet) and pion tc^ (triplet) fields defined in ( p.l| ). In principle, there 
are other dynamical fields corresponding to the remaining SU(iVj) generators (such 
as the charged pion fields n"^), but one can show that these fields do not contribute 
to the domain wall background but simply remain in their vacuum states. These 
fluctuations affect the overall energy density, but do not affect the properties of the 
domain wall such as the surface tension and so we neglect these in what follows. (See 
for example ||3^ where the Nf = 2 case is explicitly considered and vr^ = along the 



entire domain wall.) 



(4.1) 

(pT = (pu - (pd , vr° = ^^0T • 

In what follows we also need to know the masses of the relevant fields in terms of 
parameters of the effective potential (p.2|): 

< = ^, "^a = ^(l-e'), (4.2) 



TT £2 ^ a £-2 

J n J a 



and the rj' mass relation 



2 o AM , , 

= jvf/f ^ + J2- ' (4-3) 



^If one assumes that all of the rj' physics comes exclusively from the phase of the chiral condensate 
rather than from gluon condensate, then one might argue (using the linear sigma model) that the 
walls are classically unstable. Were this the case, then m,,/ ~ mq ^ in the chiral limit and 
the U(l) problem would remain unresolved 0, Thus, we see the importance of the gluon 

condensate E ^ {G'^) 3> niq which ensures the classical stability of the domain walls for iVc > 3 as 
explained in Section |^. 
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follows from (|2.15| ) with Nf = 2. Here we have neglected all possible mixing terms, 
and have introduced the following notations 



~ 2 ~ 2 

Md- Mu , , 

^ Md + M.^ ^ ^ 



4.1 Domain wall equations 

To study the structure of the domain wall we look at a simplified model where one half 
of the universe is in one ground state and the other half is in another. The fields will 
orient themselves in such a way as to minimize the energy density in space, forming 
a domain wall between the two regions. In this model, the domain walls are planar 
and we shall neglect the x and y dimensions for now. Thus, a complete description 
of the wall is given by specifying the boundary conditions and by specifying how the 
fields vary along z. 

The contribution of the light degrees of freedom to the energy density of a domain 
wall is given by the following expression^ 

a = d. + ^ + ^ + 0^, a) - Knin j (4.5) 

where the first three terms are the kinetic contribution to the energy and the last 
term is the potential. The kinetic term is actually a four divergence, but we have 
assumed the wall to be a stationary solution — hence the time derivatives vanish — 
and symmetric in the x-y plane. The only dependence remaining is the z dependence. 
Here, a dot signifies differentiation with respect to 2;: a = ^• 

Now, to find the form of the domain walls, it is convenient to use a form of the 
potential which follows from (|2.8|) : 

0T,a) = -2M (^os{^ + a^ cos^ + ^sin(^^ + a j ^) 

(4.6) 



Here we have redefined the fields 0u — 0u + a and 0^ — 0^ + a in order to remove 
the axion field from the last term ~ E and to insert it into the term ~ M. To 
minimize these equations, we can apply a standard variational principle and arrive 



^ There may be additional contributions to the wall tension from the heavy degrees of freedom 



integrated out to obtain (2.2). These will be discussed in Section 5.3 
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at the following equations of motion for the domain wall solutions: 

4^sf^ f<pT\ . f<Ps ^ \ , . (<Pt\ f^s ^ \^ E . (ct>s 

= COS — sm ha — t sm — cos h a H sm — 

8M V2 ) ^ \2 J \2 J MN, \N,^ 

(4.9) 

where the last term of Equation ([4.9|) should be understood as the lowest branch of 
the multivalued function described by Equation ( p.2|) . Namely, for (ps ^ [0,vr], this 
should be interpreted as sin((</)5 — 27rl)/Nc) with the integer / chosen to minimize 
the potential term cos((05 — 2txI)/N^. For example, with n < (ps ^ 27r, the last 
term should be of the form sin ((05 — 2n) /Nc). 

Notice the following features: first, the trigonometric terms on the right hand 
side are of, at most, order 1; thus the scale for the curvature (or rather, the second 
derivative) of the domain wall solutions is limited by fl/M and fl/M etc. In partic- 
ular, the axion domain wall must have a characteristic scale larger than m~^/(l — ^^) 
and the pion domain wall must have a scale larger than m~^. The last term in 
equation governing the (j)s field can potentially be somewhat larger than 1, hence the 
smallest scale for the ips field is related to the rj' mass. We see immediately that an 
axion domain wall must have a structure some thirteen orders of magnitude larger 
than the natural QCD scale and that the rj' field can have structure one order of 
magnitude smaller than that of the pion field. 

4.2 QCD domain walls 

Here we consider the most important case of the QCD domain wall solution which 
exists with or without an axion field. So we now set a = 0. The equations of motion 
become: 

— sm — cos — ] ^ cos — sm ' 



8M \2 J \2 J ^ \2 J \2 

hfl ( 4>t\ . { 4's\ ^ . ( 4>t\ { 4's\ , E . [ (ps 



cos( — sin( — — f sin( — cos( — H sin( — . (4.10) 

SM \2 ) \2 ) ^ \2 ) \2 ) MNc \nJ ^ ' 

For convenience, we shall label the vacuum states using the notation (0„, 0^). Thus, 
we have only one physical ground state {4>u,(pd) = (0,0), however, because of the 
conserved topological current (|3.1| ), classically stable domain walls can form and in- 
terpolate from the ground state (0„, 0^) = (0, 0) along a path which is not homotopic 
to the null path. To classify the paths we use the redundant notation where (0, 0) 
and {2tt, 0) etc. are considered as different states and we talk about the field inter- 
polating between these states. Keep in mind that this is only a way of classifying 
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the homotopy classes and that in fact all the states represented by (27m, 2iTm) for 
integers m and n are one in the same vacuum state. 

The simplest domain wall is described by a continuous transition from the ground 
state {(pu, (pd) = (0, 0) to the state labelled 0^) = (2vr, 0) as described by the 
vacuum solution Equation ( |2.11| ) with 6 = 27i (or, equivalently, Z = — 1 in Equa- 
tion ( |2.2| )). This wall corresponds to a single winding around the U(l) manifold. It 
is also possible to wind in the opposite sense. To summarize, the two topologically 
stable domain walls of minimal energy correspond to one winding in each direction 
and are classified by the transitions from (</>„, cpd) = (0, 0) to: 

Soliton (0„,0d) = (27r,0). 
Antisoliton {(pu^'t'd) = (— 27r,0). 

The solutions which wind from (0, 0) to (0, ±27r) are not topologically distinct from 
these due to the other pion fields and have a higher energy. Note, however, that in the 
chiral limit m„ = and ^ = 0, thus the transitions to ((/)„, (pd) = (0, ±27r) have the 
same energy and there is a degeneracy. If rriu > rrid, then these transition in (pd are 
the minimal energy solutions and the (pu solutions above become unstable. In reality 
rrid > rriu, and the transitions to {(pu, (pd) = (=t27r, 0) are the only stable transitions. 

The general case of Equations ( [4.10|) cannot be solved analytically and we present 
the numerical solution of Equation ( fl.lOp in figure 0. In order to gain an intuitive 
understanding of this wall, we examine the solution in the limit rriTr -C rrirj'- In this 
case, the last term of ( |4.1CI| ) dominates unless (ps is very close to the vacuum states. 



~ 27m. Thus, the central structure of the (ps field is governed by the differential 



equation 



8E 



Now, there is the issue of the cusp singularity when (ps = n because we change 
from one branch of the potential to another as expressed in Equation ( p.2| ). By 
definition, we keep the lowest energy branch, such that the right hand side of 
Equation ( |4.11| ) is understood to be the function sm{(ps/Nc) for < 05 < vr and 
sin((05 — 27r)/A'^c) for 7i < (ps < 27i. However, we notice that the equations of motion 
are symmetric with respect to the centre of the wall (which we take as z = zq), hence 
(ps = only at the centre of the wall and not before, so we can simply look at half of 
the domain, z G (— c>o, 0], with boundary conditions (ps = Q sX z = — oo and (ps = it 
at z = 0. The rest of the solution will be symmetric with (ps = 27t at z = +oo. 
Equation ( [4.11| ) with the boundary conditions above has the solution (recall that 



= (pu + 

' ANctan-^ 



t>s{z) =(pu + 



tan ^ e'^^^-^") 



2tt - ANr tan 



-1 



tan ^ e'^^^^-^o) 



z <Zo 

(4.12) 

Z>Zo 
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where Zq is the position of the centre of the domain wall and 

2y/2E 

/i = — — — , lim u = m,^> (4.13) 

is the inverse width of the wall, which is equal to the m^/ mass in the chiral limit 
(see Equation ( [4. 3D ). Thus, we see that the dynamics of the central portion of QCD 
domain walls is governed by the rj' field. We shall also refer to the (ps transition ( [4.12| ) 
which occurs in several places (see for example Section [4.3.2| ) as the rj' domain wall. 



The first derivative of the solution is continuous at z = zq, but the second derivative 
exhibits a finite jump. 

Before we continue our discussions regarding the structure of QCD domain walls, 
a short remark is in order: The sine-Gordon equation, which is similar to Equa- 
tion ([4.1 1|) with a cusp singularity, was first considered in where a solution similar 



to ( |4.12| ) was presented. There is a fundamental difference, however, between domain 



walls discussed in and the domain walls we consider here. In [Q, the domain 
walls were constructed as auxiliary objects in order to describe a decay of metastable 
vacuum states which may exist under the certain circumstances. The walls we dis- 
cuss here classically stable physical objects where the solutions interpolate between 
the same vacuum state; their existence is a consequence of the topology of the U(l) 
singlet T]' field. This topology, represented by the exact symmetry 9 = 6 -\- 27in in 
Equation ( p.2|) is a very general property of QCD and does not depend on the specific 
choice of parameters or functional form of the effective potential. Similar equations 
and solutions with application to the axion physics were also discussed in [^ . 

The solution described above dominates on scales where \z\ < however, 
the isotopical triplet pion transition can only have a structure on scales larger than 
^ and so the central structure of the t]' wall can have little effect on the 
pion field. Indeed, we can see that, for \z\ ^ /x^^, 05 is approximately constant with 
the vacuum values. Making this approximation, we see that the isotopical triplet 
field is governed by the equation 

,2 „■ / (PT 



2m^sin(:^). (4.14) 



This has the same form as ( [4.11|) and hence the solution is (recall that 0t 



fstan-^Ttanf e"'-(^-^")l , z <t: Zq - U~\ 

[27r-8tan-i[tanf e-"^-(^-^o)] , z Zq + fi~\ 

which is a reasonable approximation for all z. Numerical solutions for the (ps and 0t 
fields are shown along with the same solution in terms of the and (pa in figure ^ 
As we can see from the explicit form of the presented solution, the rj' transition is 
sandwiched in the pion transition. This is a key feature for some applications of this 
type of the domain wall as discussed in H] for example. 
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Figure 3: Basic form of tlie QCD domain walls. Notice that the scale for the pion 
transition is larger than for the r]' transition and that the width of the r/' wall is set by the 
scale rrir)'- 



The contribution to the wall surface tension defined by Equation (|4.5| ) can be 
easily calculated analytically in the chiral limit when the analytical solution is known 
and is given by Equations ( |4.12|) and ([4.15|) . Simple calculations leads to the following 
contribution from the pion and rj' fields. 



In case when rriq 7^ 0, an analytical solution is not known, but numerically, a is close 
to the estimate (|4.16|) . 

4.3 Axion dominated domain walls 

In the previous subsection when the QCD domain walls were discussed, the axion was 
not introduced as a dynamical field. In this subsection we assume that the axions 
exist. In this case there are domain walls in which the axion is the dominant player. 
The introduction of axions, in most cases, makes the domain wall an absolutely stable 
object. Our case is no exception and the axion model under discussion, (which is 
the N = 2 axion models according to the classification is an absolutely 

stable object. At the same time it is well-known [^, that stable domain walls 
can be a cosmological disaster. We do not address in this paper the problem of 
avoiding a domain wall dominated universe. Rather, we would like to describe some 
new elements in the structure of axion domain walls, which were not previously 
discussed. 

The first and most natural type of the axion domain wall was discussed by Huang 



and Sikivie p^] who neglected the rj' field in their construction. We shall refer to 



this wall as the Axion-Pion domain wall (a-^). As shown in 1391, it has a width of 
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the scale ~ m^^ ^ ^qcd both the axion and vr meson components. As Huang 
and Sikivie expected, the t]' plays a very small role in this wall. In what follows we 
include a discussion of this type of domain wall for the completeness. 



Our original result is to de- 
scription a new type of the axion 
domain wall in which the t]' field 
is a dominant player. We shall 
call this new solution the Axion- 
Eta' domain wall (a^/). This new 
type of the domain wall was con- 
sidered for the first time in 

possible source for galactic 
magnetic fields in early universe. 
In what follows we give a detail 
description of the solution for the 
a^/ wall. Here we want to men- 
tion the fundamental difference 
between the wall discussed in 
and the a^/ wall introduced in ||^ 





-5 



a(z) 





zma 



igure 4: Basic form of the a-ir domain wall. Notice 
hat all the fields have a structure on the scale of 
and that the isotopical singlet fields plays a very 
small role. 



Unlike the a^^ wall which has 
structure only on the huge scale 

of m~^, the ttrf' wall has nontrivial structure at both the axion scale as well as 
at the QCD scale m~,^ ~ ^qcd- '^^^^ reason for this is that, in the presence of the 
non-zero axion field (which is equivalent to a non-zero 6 parameter), the pion mass 
is efficiently suppressed due to its Goldstone nature, thus the pion field follows the 
axion field and has a structure on the same scale. The 1]', however, is not very 
sensitive to 6 and so it remains massive. 

Again, the a^/ solution has a 
sandwich structure with the sin- 
glet transition occurring at the centre of the wall. One can adopt the viewpoint that 
the Qri' domain wall is an axion domain wall with a QCD domain wall sandwiched 
in the centre. This phenomenon is critical for applications involving the interac- 
tion of domain walls with strongly interacting particles. Indeed, there is no way 
for the wall to trap any strongly interacting particles, like nucleons, because of 
the huge difference in scales; the a,,/ wall, however, has a QCD structure and can 
therefore efficiently interact with nucleons on the QCD scale Aq^j-,. Regarding the 
wall tension ( [4.5| ), it is dominated by the axion physics and thus has the same or- 
der of magnitude for both types of the axion domain walls which is proportional to 
a ~ M/nia ~ fafnfnn as found in 



20 



4.3.1 Axion-pion domain wall 

The solution discussed by Huang and Sikivie corresponds to the transition (a, 0^) : 
(0, 0, 0) — >■ (tt, tt, — tt), i.e., by a transition in the axion and pion fields only. This tran- 
sition describes the a^^ wall. Indeed, in terms of (a, 0^, 0^) this transition corresponds 
to a nontrivial behaviour of the axion and triplet pion component 0^ of the 0^1, 0^ 
fields: (a,05,0T) : (0,0,0) — > (vr, 0,27r). The other transition (the a,,' wall) corre- 
sponds to the transition (a, 0^, 0^) : (0, 0, 0) (— tt, tt, vr) in which the singlet rj' 
field is dominant: (a, 05, 07-) : (0, 0, 0) (— vr, 27r, 0). It will be considered later on. 

Huang and Sikivie discussed the solution to this wall in the limit where the rj' 
field is extremely massive and hence they neglected its role. It can be integrated out 
which effectively corresponds to fixing it ri'{z) = 0. Indeed, if we simply fix ri'{z) = 
in our equations, then we reproduce their solution. When E is large, then m^r <C rrirj' 
as Huang and Sikivie assumed, the effects of the rj' particle can be neglected and the 
solution for the the axion and pion fields presented in is valid for the boundary 



conditions described above. We plot this numerical solution which includes the t]' 
effects in figure ^ As announced above, the solution has the only scale of ~ 
for both components, axion as well as the pion field ~ 0t. The i]' field remains close 
to the its vacuum value and only slightly corrects the solution. 

4.3.2 Axion-eta' domain wall 

Having looked at the solution of the a^^ wall, we now investigate the structure of the 
a^/ domain wall which is a new solution. This solution corresponds to the transition 
(a, 0„,0rf) : (0,0,0) (7r,7r,7r). Now the singlet rj' field undergoes a transition 
instead of the triplet pion field: (a, 05, 0t) : (0, 0, 0) — > (— vr, 27r, 0). As we discussed 
above, the singlet field never becomes massless, and therefore, a new structure at the 
QCD scale ~ ^qcb emerges in sharp contrast to the well-studied a^, domain wall [p9 



where no such structure appears. We should note that the potential ([4.6|) has the 
same vacuum energy at 05 = as well as at 05 = 27r due to the change of 05 field 
to the lowest energy branch at 05 = vr as described by Equation ( p.2|) . Therefore, 
this domain wall interpolates between two degenerate states, and thus, like the a^r 
domain wall [^, the solution under these considerations is absolutely stable. 



For \z\ ^ the last term in ( [4.9| ) is negligible and so the solution behaves like 
the wall. What happens is that, away from the wall, the axion field dominates and 
shapes the wall as it does in with the solution. Again, the pion mass is suppressed 
and 05 ^ 0. As z ~ fi~^ however, the last term of (|4.10| ) starts to dominate the 
behaviour. At this point, the a^/ wall undergoes a sharp transition similar to the 
QCD domain wall described by ([4.12|) . We plot this solution along with a blowup in 
figure 1^. Notice also, that the singlet field cancels the effects of the axion near the 
center of the wall, and so the pion field becomes massive again as it undergoes its 
transition. 
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Figure 5: Basic form of tlie a^' domain walls with a closeup where the axion field a ~ 7r/2. 
Notice that the large scale structure is similar to that of the wall, but that there is also 
a small scale structure on the scale of m^/. Near the centre of the wall, the pion regains its 
mass and undergoes a transition on the scale of m.,^. 



5. Decay of the QCD domain walls 

We do not have much new to say regarding the generalities of domain walls, nor do 
we have a resolution for the general problem of avoiding a domain wall dominated 
universe. We have nothing new to say regarding the stability or evolution of axion 



domain walls |^-|^ (which were also discussed in the previous sections), see ^0 



and references therein for a recent review on the subject. We refer the reader to 



the nice text book |Q for general discussion about domain walls, other topological 
defects and their role in the early universe. This section is devoted specifically to the 
QCD domain walls discussed in Section ^72. 



Up to now, we have treated the domain walls as topologically stable objects. If we 
only consider the low energy degrees of freedom present in the effective theory (|2.2| ) 
then this is certainly the case as discussed in Section We shall also argue that this 
picture of classically stable domain walls is correct in the appropriate large Nc = oo 
and chiral = limit. In reality, however, one must consider the heavy degrees of 
freedom that were integrated out to obtain ( ^.2[ ). For finite Nc, one finds that QCD 
domain walls (not axion domain walls) are unstable on the quantum level due to a 



tunnelling mode of decay as described in Section |3.1| . We estimate the relevant life- 
time of these walls with respect to this quantum transition and show that, although 
the walls don't live long on a cosmological scale, exponential suppression of the 
tunnelling decay mode may result in the lifetime of the walls being much larger than 
the QCD scale Aq^j-,. We argue that this exponential suppression remains in the 
large Nc limit where we have theoretical control, and may also occur numerically for 
physical values of the parameters. Therefore, these walls do not pose a cosmological 
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problem as one might naively suspect/ 

Even though QCD domain walls are ultimately unstable, they may still play 
an important role in physical processes with timescales comparable to the lifetime 
of the walls. QCD domain walls are not likely to exist today, however, they may 
play an important role shortly after the QCD phase transition in the evolution of 
the early universe, or in heavy ion collisions after the transition from a quark-gluon 
plasma to the hadronic phase when system cools. In these cases, QCD domain walls 
with a finite lifetime may play an important role. Indeed, according to the standard 
theory of cosmological phase transitions [0, |3^ , if, below a critical temperature 



the potential develops a number of degenerate minima, then the choice of minima 
will depend on random fluctuations in fields. The minima that the fields settle to 
can be expected to differ in various regions space. If neighbouring volumes fall into 
different minima, then a kink (domain wall) will form as a boundary between them.^ 
The relevant question becomes: "Is the lifetime of QCD domain walls large enough 
to be of physical interest?" In what follows, we hope to convince the reader that the 
answer to this question may in fact be: "Yes!" . 

5.1 Estimating the decay rate 

In the following, we estimate the lifetime of the QCD domain walls due to a tunnelling 
process whereby a hole forms in the domain wall and expands, consuming the wall. 
A quantitative calculation of this rate is presently beyond our control, but we can 
estimate the magnitudes of the effect through a semi-classical approximation. In 



Section 2^ we shall argue that these approximations are valid asymptotically in the 
large Nc hmit and that the decay rate is exponentially suppressed. 

The decay mechanism is due to a tunnelling process which creates a hole in 
the domain wall which connects the (pd) = (0,0) domain on one side of the 
wall to the (</>«, 0d) = (2vr, 0) domain on the other (see Equation ( p.llp ). Passing 



through the hole, the fields remain in the ground state. This lowers the energy of the 
configuration over that where the hole was filled by the domain wall transition by an 



^We should remind the reader once again that the existence of the QCD domain walls described 
above is a consequence of the well-understood symmetry 6—^6 + 27m and is a consequence of the 
topological charge quantization; their existence is not based on any model-dependent assumptions 
we have made to support the specific calculations in the previous sections. The question now is not 
the existence of these walls, but whether or not they live long enough to affect relevant physics. 

* Whether or not domain wall actually will form at the QCD transition in the early universe is 
an unresolved issue. Presently it is believed that the QCD phase transition is actually a smooth 
crossover (see for a review of the QCD phase diagram). If this is the case, then it is possible 
that no domain walls form because the universe cools very slowly. To resolve this question, one 
must estimate the relaxation timescales involved: if the crossover is sharp enough, then domain 
walls may still form. At the RHIC, however, the system is quenched due to the rapid expansion of 
the colliding ions and so the formation of domain walls is very likely if an appropriate part of the 
phase diagram is explored by the reaction. 
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amount proportional to where R is the radius of the hole. The hole, however, must 
be surrounded by a string-like field configuration. This string represents an excitation 
in the heavy degrees of freedom and thus costs energy, however, this energy scales 
linearly as R. Thus, if a large enough hole can form, then it will be stable and the 
hole will expand and consume the wall. 

This process is commonly called quantum nucleation and is similar to the decay 
of a metastable wall bounded by strings, and we use a similar technique to estimate 
the tunnelling probability. The idea of the calculation was suggested by Kibble 



and has been used many times since then (see the textbook |38] for a review). The 
most well known example of such a calculation is the calculations of the decay rate 
in the so-called Npq = 1 axion model where the axion domain wall become unstable 
for a similar reason due to the presence of axion strings [^-[Q. However, as 
was emphasized in [^, the existence of strings as the solutions to the classical 



equations of motion is not essential for this decay mechanism (see below). Some 
configurations, not necessarily the solutions of classical equations of motion, which 
satisfy appropriate boundary conditions, may play the role played by strings in the 
NpQ = 1 axion model. 

To be more specific, let us consider a closed path starting in the first domain 
with (0„, 0(i) = (0,0), which goes through the hole and finally returns back to the 
starting point by crossing the wall somewhere far away from the hole. The phase 
change along the path is clearly equals to (pu + 4>d = 27r. Therefore, the absolute 
value of a field which gives the mass to the 1]' field (the dominant part of the domain 
wall) has to vanish at some point inside the region encircled by the path. By moving 
the path around the hole continuously, one can convince oneself that there is a loop 
of a string-like configuration (where the absolute value of a relevant field vanishes 
such that the rj' singlet phase is a well defined) enclosing the hole somewhere. In 
this consideration we did not assume that a hole, or string enclosing the hole, are 
solutions of the equations of motion.^ They do not have to be solutions. 

However, if we want to describe the hole nucleation semi-classically |53|, ^ 



then we should look for a corresponding instanton which is a solution of Euclidean 



^It is quite obvious that such a configuration cannot be described within our non-hnear a model 
given by Equation ( |2.2|) where it was assumed that the gluon as well as the chiral condensates are 
non-zero constants. In this case, the singlet phase is not well defined everywhere. However, in the 
case of a triplet tt meson string, such a configuration can easily be constructed within a linear cr 
model by allowing the absolute value of the chiral condensate to fluctuate along with the Goldstone 
phase (tt meson field). The a term in the linear a model essentially describes the rigidity of the 
potential. Indeed, the corresponding calculations within a linear a model were carried out in ||5^ 
where it was demonstrated that the solution describing the tt meson string exists, albeit unstable 
as expected from the topological arguments. To carry out a similar calculations in our case for 
the singlet rj' phase, one should allow fluctuations of the gluon fields: the fields that give mass to 
T]' meson and that describe the rigidity of the relevant potential. (See ( |3.1[ ) and the surrounding 
discussion.) 
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(imaginary time, t = it) field equations, approaching the unperturbed wall solution 
at r — s> ±00. In this case the probability P of creating a hole with radius R per area 



S per time T can be estimated as follows |5^, |5^, |57 

3 



P 
'ST 




e-^oxDet, (5.1) 



where 5*0 is the classical instanton action; Det can be calculated by analyzing small 
perturbations (non-zero modes contribution) about the instanton, (see for an 



explanation of the meaning of this term) and will be estimated using dimensional 
arguments; and (\/5'o/(27r)) is the contribution due to three zero modes describing 
the instanton position. 

If the radius of the nucleating hole is much greater than the wall thickness, we 
can use the thin-string and thin-wall approximation. (The critical radius Rc will be 
estimated later and this approximation justified). In this case, the action for the 
string and for the wall are proportional to the corresponding worldsheet areas 



So = A'KR^a - -TiR^a . (5.2) 
3 

The first term is the energy cost of forming a string: a is the string tension and AttR"^ 
is its worldsheet area. The second term is energy gain by the hole over the domain 
wall: a is the wall tension and A/SirR^ is its worldsheet volume. The world sheet 
of a static wall lying in the x-y plane is the three-dimensional hyperplane z = 0. In 
the instanton solution, this hyperplane has a "hole" which is bounded by the closed 
worldsheet of the string. 

Minimizing (|5.2| ) with respect to R we find the critical radius 

S„ = l^. (5.3) 

The lorentzian evolution of the hole after nucleation can be found by making the 
inverse replacement r — > —it from Euclidean to Minkowski space-time. The hole 
expands with time as x"^ + y"^ = R^ + t"^ , rapidly approaching the speed of light. 

To estimate the appropriate string and wall tensions, we must step back from the 
effective theory (|2.2|) and include the heavy degrees of freedom that allow the fields 
to tunnel. Unfortunately, to provide a well justified estimate of these parameters in 
standard QCD is very difficult because we do not know how to quantitatively include 
the effects of the heavy degrees of freedom. We shall argue, however, that we regain 



-'^^The estimate given below is designed for illustrative purposes only, and should be considered 
as a very rough estimation of the effect to an accuracy not better than the order of magnitude. 

^^The three zero modes in our case should be compared with the four zero modes from the 
calculations of This difference is due to the fact that in the decay of three dimensional 
metastable vacuum state was discussed. In our case, we discuss a decay of a two-dimensional object. 
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this theoretical control in the chiral and large Nc limits. We postpone this discussion 
until Section |5.2| , but we summarize the Nc dependence of these parameters here: 



Nc<a, VNc<a. (5.4) 

These lower bounds for the string tension a and wall tension a are well justified. 
For our argument, however, we need an upper bound on the wall tension. As we 
shall argue later, we expect that the lower bound for a in ( ^.41) is actually the upper 
bound and that: 

(j~x/iVc- (5-5) 

However, there is some speculation about the correct upper bound and one might 
argue that that the upper bound could be of order A^^. While we strongly suspect 
that (|5.5| ) is correct, we cannot prove this and thus consider the possibility of larger a. 

In either case, from (|5.3| ) we see that, in the large limit, the probability of 
producing a hole (|0| ) is exponentially suppressed by the factor of at least e""^" ~ 
g-A^c jj^ ^Y^Q (.g^gg of 0" ~ A'^c (In this case, our semiclassical estimate of 5*0 is not 
numerically justified, however, we believe that it is parametrically still valid as we 
shall discuss in Section p75|) . We believe, however, that ( p.5| ) holds and thus that the 
actual suppression is much stronger 

e-^o ~ e"^' . (5.6) 

The only way to kill the exponential suppression is to arrange for a > Nc which, as 
we shall discuss, has no phenomeno logical or theoretical support. Thus, in the limit 
A'c — >■ oo, we believe that no tunnelling is supported and the domain walls become 



stable. In Section |5.6| we shall extrapolate these results to the A"c = 3 limit and, 
although we lose the theoretical control gained in the large Nc limit, we argue that 
the qualitative picture might remain the same. 

5.2 Heavy degrees of freedom 



As discussed in Section ^TT], we cannot properly discuss tunnelling with Equation ( |2.2| ) 
where the gluon degrees of freedom are integrated out and replaced by their vacuum 
expectation values: such a theory cannot describe the strings which are responsible 



for the domain wall decay described in Section [Ol . Instead, we must take one step 
back and describe the dynamics of the gluon condensate by considering the original 
Lagrangian |]12|, |13| which includes the complex gluon degrees of freedom h. In this 
effective theory, the complex gluonic field will sit in a Mexican-hat potential with 
\{h)\ = const. The string represents a localized region about which h = 0. In the 
effective theory ( p.2|) this represents a singularity, but here the singularity is allowed, 
although it is energetically costly. The string tension is associated with the energy 
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cost required to form the string and this is directly related to the height of the peak 
of the Mexican-hat potential. 



The effective potential of the original Lagrangian [|1^, |T3| which includes the 
complex gluonic degrees of freedom h is given by:^^ 



v{e,h,u) 



h 

AN, 



log 



h 



2eE 



DetU 



- iTiMU 1 +H.C. 



(5.7) 



This satisfies all the conformal and chiral anomalous Ward identities, has the correct 
large Nc behaviour etc. (see [0, |13|, |3^ for details). Integrating out the heavy h field 



will bring us back to Equation (|2.2| ). (It is interesting to note that the structure of 
Equation ( ^.7|) is quite similar in structure to the effective potential for SQCD 
and gluodynamics p8| .) 

Let us stop for a moment to consider the form of this potential. When we 
integrate out the heavy gluon degrees of freedom, we replace h with its vacuum 
expectation value {h). The vacuum expectation value of the h is given by 



such that 



{h) = 2Eexp{-i^(f)i/N^) = 2Ee~'^'^^' 
becomes 



V{U) 



h={h) 



-E cos 



cos ( 



(5. 



(5.9) 



in agreement with ( p.8|) . With h fixed at its vacuum expectation value, the singlet 
combination (ps = exhibits the U(l) topology described above and our domain 
walls are stable. Now, however, we allow the gluon condensate to fluctuate. We 
parameterize these fluctuations in polar coordinates by a radial component p, and 
an angular component (: 



h = pe''^^{h) = pexp 



4 ' ■ 



(5.10) 



Here is the decay constant and 0^ is the dimensionless phase angle. The fields p 
and ( are both heavy, real, physical fields. In the large Nc limit, their masses are 
much heavier than the pion or t]' masses: ~ rup 3> m^/ . The role of the ( field will 
be important in explaining the physics of the cusp singularities in ( |2.2|) and will be 
discussed in Section |5.4| . For now we set 6 and 0^ to zero. In terms of the remaining 
physical fields, the potential V becomes 



V{p, <Ps) = +Ep{\og p-l) cos 



9s_ 
N. 



(5.11) 



^^Note: the expression ( 5.7) is valid only for one branch. For a more accurate expression and treat- 
ment of the subsequent minimisation process which carefully accounts for the different branches, 
see the original paper Isl] . 
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where we have neglected the terms proportional to M since they only contribute a 
constant offset due to the 0t field. For an early phenomenological discussion of the 
potential ( |5.11| ) without the (j)s fields, see . 



Now the combined degrees of freedom p and (ps are no longer restricted to the 
circle p = 1 as they were in the effective theory ( |2.2| ). The U(l) topology is no longer 
a constraint of the fields and thus, the walls are not topologically stable. Instead, 
the restriction of p ~ 1 is dynamical, and made by a barrier at h = 0. Thus, with 
these degrees of freedom, the fields parameterize the plane, however, the potential 
is that of a tilted Mexican-hat with a barrier at /i = 0. The barrier is high enough 
that a domain wall interpolating around the trough of the hat is classically stable. 
If the barrier were infinitely high as we assumed when we fixed p = 1 as we did 
in ( p. 21) , then the (ps field could wind around the barrier and would be topologically 
stable. With a finite barrier, however, the field can tunnel through the barrier as 
described above. This situation is analogous to the case of the string and peg shown 
in figure 0. We show a more accurate picture^'^ of the barrier ( |5.11| ) in figure |^. The 
relative heights of the the peak and troughs are given below: 

AV^Peak = AVi + A^s = E ~ iV^ (5.12) 
AVi = Ecos ^ ~ , (5.13) 

Al^s = ^ ( 1 -cos^ ) ~ 1. (5.14) 



We have emphasised the Nc dependence here because the effective theories are really 
only well justified in the large Nc limit. 

The most important property of the potential is the following: The absolute 
minimum of the potential in the chiral limit corresponds to the value V^i^ = —E 
which is the ground state of our world with p = 1 and 0^ = 0. At the same 
time, the maximum of the potential (|2.2|) , where one branch changes to another 
one is V = —E cos{7t /Nc) where 0s = (we are still taking 6 = 0). This corre- 
sponds to the point p = 1 and 0s = vr in the potential (|5.11|) . Thus, the trough 



of the Mexican-hat is given where h = (h), i.e. at radius p = 1 and the maxi- 
mum V = —Ecos{tt/Nc) of the potential ( p.2| ) is exactly the barrier through which 
the rj' field interpolates to form the QCD domain wall. It is important to note 
that the height of the barrier for the potential ( |2.2| ) is numerically is quite high 
~ E{1 — cos-^), but vanishes in large Nc limit. Indeed, in this limit, the peak of 
the barrier is degenerate with the absolute minimum Vram = —E of the potential as 
it should be.^'^ The height of this barrier describes how much the Mexican-hat is 
tilted. 



^^We are indebted to Misha Stephanov for suggesting this nice intuitive picture for explaining 
the domain wall decay mechanism. 

^^Remember, the if direction becomes flat in the large Nc limit as m,y 0. 
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AVo ~ 1 




Figure 6: Profile of the "Mexican-hat" potential ( p.ll| ). The slice is made along the axis 
through 05 = vr to the left and (/^^ = to the right. The trough of the potential lowers from 
the cusp at p = 1, = vr where V = —E cos{tt /Nc) down to the vacuum state Vmin = —E. 
The hump where y = is at the origin is where h = pexp{i(ps/Nc) = and hence the 
singlet field (ps can have any value at this point. It is by passing across this point that a 
QCD domain wall can tunnel and a hole can form. The important qualitative features of 
this potential are the height of the central peak AVpcak ( |5.12| ), the relative heights of the 
cusp to the peak AVi (5.13) and the trough to the cusp AV2 (|5.14|). 



The other important property of the potential ( |5.11 ) is its value where the singlet 
phase 'pi is not well defined at the peak of the Mexican-hat. From ( |5.11| ) it is clear 
that this occurs for /i = 0. It is this peak of the Mexican-hat in figure |^ (or the peg 
in figure ^ that classically prevents p ^ and that makes the QCD domain wall 
classically stable. When a hole tunnels through the wall, it must be surrounded by 
a "string" where the field passes through the region h = 0. The height of the peak 
thus contributes to the energetic cost of creating such a string (and hence the hole). 
The potential ( |5.11| ) vanishes at the centre of the string V{h = 0) = 0, which implies 
that the barrier at /i = is quite high: AVpeak = E. As expected, the barrier at 
h = should be order of ~ in contrast with the barrier to the rj' domain wall 
where one expects a suppression by some power of Nc- We also note that the total 
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number of distinct classically stable solutions can be estimated from the condition 
that cos{7rk /Nc) > where the barrier for the rj' field is still lower than the peak 
h = 0. Thus, —Nc/2 < k < Nc/2 where k labels the winding number of the solution. 
Thus we see that for Nc <2 there is no admissible classically stable solution, but for 
Nc > 3 classically stable solutions are allowed. In our model with potential 
we see explicitly that, for Nc = 3, there is one classically stable domain wall while 
for Nc = 2 there are no classically stable solutions. 



5.3 String and wall tensions 



In order to make the semi-classical estimates of the decay rate (see Equations (|5.1j ) 
and (|5.3|) ) we must estimate the string tension a and the wall tension cr. We start 
with the string tension a. 

Within the effective theory ( p.7|) true strings are not supported because the U(l) 
symmetry is broken by the anomaly. In the large Nc limit, however, this symmetry 
is restored (and the 1]' becomes massless). Thus, it makes sense to consider global 
U(l) strings in the large Nc limit. To properly estimate a at finite Nc one must 
numerically minimizes the energy of a configuration with a domain wall bounded by 
a string-like configuration. This could be done numerically, but is beyond the scope 
of the present paper. In the large Nc limit, the estimates presented here become 
reliable. 

To estimate the string tension a, we first consider an isolated global string in a 
flat Mexican-hat potential 

£=^9^<l>*9^'$-y(|$|), (5.15) 

where $ is a generic complex scalar field that will be identified with the glue-ball 
field h discussed in Section p.2\ . To match l^(|$|) with ( [5. 7] ) we must take the limit 
Nc ^ oo so that the U(l) symmetry $ e*°<l> (equivalently h e*"/i) is restored. 
In this string lying along the 2;-axis with winding number n will be described 

by the complex field configuration 

$(r, 6) = p{r, e)e^'^("'^) = p(r)e*"^ = p(r) exp (^in tan'^ ^) (5.16) 

in the x-y plane. The exact radial dependence p{r) will be such as to minimize the 
energy density along the string and will minimize the energy density 



a 



j I {]fif3iP + ^p%0,^ + V(p) - J dxdy , (5.17) 
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We have assumed here that the fields p and are canonically normahzed. 

Within some radius r^, the radial dependence of the string will vary from to 
(p) and there will be a core energy contribution to the string tension 

«o = vrr,Vo, (5.19) 

where Vq represents an average core energy. Far away from the string, p will assume 
its vacuum expectation value and the total string tension (energy density per unit 
length) will be 

a~ao + 7r/ — ^dr ~ Ocore + 7rn^(p)^log-R, (5.20) 



r 



where we have absorbed all contributions from the region of size into the constant 
Ocore- For an isolated string, the string tension is infinite, but in most cases, the lateral 
extent of the string is limited by some upper radius R that must be determined by 
the dynamics of the strings. In our case, the string is embedded in the end of a 
domain wall, so the relevant scale for R will be the wall thickness. 

Far from the core of the string, the only degree of freedom is the dimensionless 
phase 0: 



$ = (p)e^'^. (5.21) 



Thus, the kinetic term is 



]^{pfd,cpdy. (5.22) 



The normalization of the field p is such that this term reduces to the appropri- 
ate canonical kinetic term for the appropriate Goldstone field described by Equa- 
tion ( ^.15|) . In our case, the string is composed of the rj' field and the relevant phase 
(f) interpolates from to 27r connecting the two sides of the QCD domain wall. For 
this 7]' string, the the relevant phase is = 2ri' / (frji \/Nf) is given in ( p.l|) : thus 
(p)^ = Nff'^ijA, and the tension is: 

a ~ "core + -Nff'^iv? logi? ~ a ~ Ocore + -j^fl log R , (5.23) 

where we have set Nf = 2 and n = 1 (single winding) in the last equation, as we 
have been doing. 

To determine Ocore, one must actually minimize the string tension. This requires 
full knowledge of the potential V{p) and the nature of the field p. We do not have 



15 



To justify the use of ( [5.201 ) the radius of the string must be much smaller than the wall 



thickness. We shall show that this is the case in (5.25) 
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this information. In general, however, we expect that the scale of the core is set by 
the mass of the heavy field p 



~ m;^ ~ 1 , (5.24) 



and thus, in the large Nc limit, the string core size becomes much smaller than the 
domain wall thickness (|4.13|) 

r, ~ m"^ ~ K /i"^ ~ m"/ ~ 0Vc, (5-25) 



so one can think of the string as embedded in the edge of the domain wall. In this 
environment, the outer radius of the string R is the same order as the thickness of 
the domain wall. This justifies the approximation of the string tension ( [5.23| ) as that 
of a free string with outer radius R: if the string core had been of comparable size to 
the domain wall, then this approximation would be inaccurate, and we would have 
had to minimize the energy of the combined system of a domain wall bounded by a 
string. This situation is much more difficult to solve because the "string" would no 
longer have cylindrical symmetry. 

Unlike the tension of a global string, which has a logarithmic contribution from 
the bulk ( ^.23| ), the domain wall tension (c.f. Equation ( |4.5| ) come exclusively from 
the core dynamics. Due to the quadratic kinetic terms, there is an equipartition of 
energy and the tension is essentially 

a ~ A^K) , (5.26) 

where /S.z is the wall thickness and Vq ~ Vcove — Knin is the average potential en- 
ergy near the core. For most domain walls, the thickness is governed by the mass 
of the relevant field. Thus, for QCD domain wall tension, we have the following 
contributions from the pion and rj' transitions respectively (see potential ( |4.6| )): 

a, ~ m-'2M ~ m;'m,|(*>I'>| ~ f;m^ , (5.27) 
a,,~V^^(l-coB^)~^, (5.28) 



which agree qualitatively with the precise calculations ( 4.16| ). Thus, in the chiral 
limit rUg 0, only the rj' contribution remains relevant. Numerically, even with 
finite quark masses, the t]' contribution is much larger than the pion contribution. 
There is another contribution, however, due to heavy degrees of freedom which we 
estimate in the next section. 

5.4 Cusps 

As was emphasized in ||6^ for supersymmetric QCD, the presence of a cusp singularity 



in the effective potential ( p.2| ) can indicate that heavy degrees of freedom are playing 



32 



a role in the physics. As far as the low energy degrees of freedom are concerned, (the 
pions and rj'), the effective potential ( |2.2| ) is valid on the scales associated with these 
degrees of freedom. If we were to properly include the heavy degrees of freedom, we 
would find that the cusps are actually smooth, but only on scales comparable to the 
mass of the heavy particles rrih ^ m^r. In the centre of the domain wall, where the 
potential has the cusp, what is really happening is that the heavy fields are making 



a rapid transition with a scale length of Az ~ 



1 



Analogously to [Q, the contribution of this heavy transition to the domain wall 
tension can be estimated by (|5.26|) using the heavy fields introduced in (|5.1CI| ). The 



mass scale sets the width for the transition, but to determine the scale for the 
potential Vq requires more work. One can estimate this from the mass term of the ( 
field as follows. 

In order to correctly allow the rj' field to interpolate from 6 = to 6 = 27t, the 
field ( enters the picture as a phase in the following combination with the t]' and 6 
angle: 

exp I— + i—= h I— = exp I — ^ + i—= . (5.29) 



In order to allow the rj' field to switch from one branch of the potential to another, 
( must shift in such a way that 6^6 + 27r. Thus, with our definition, ( makes the 
transition from to 27r/^/A^c- Since there is an equipartition between the kinetic and 
potential energies, we can estimate the contribution to the wall tension from the ( 
field by using the kinetic term: 

The numerical contribution of this term to the wall tension might be quite large 
since it is proportional to the heavy mass m^. Expanding ( ^.29| ) and noting that 
the potential will depend on some energy scale ~ related to gluonic physics, we 
estimate the mass term as 

ml ~ Ar2//2 . (5.31) 

Using the standard assumption that masses do not depend on Nc, we have that 
~ 1. Thus, ~ Nc and 

(T^ ~ 1 . (5.32) 

Thus, the presence of cusps in an effective theory signals that heavy degrees of 
freedom might play an important quantitative role through contributions like a^, in 



agreement with the conclusion reached in |60|. If we make the same estimates for 



the 7]' meason, we get a.^/ ~ yiVc instead of ( |5.32| ) due to the unique way that the 
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Vj' couples to the Q parameter through the term {Q — i logDet \J\ and the unique 
dependence: m^/ ~ iVc ^^"^ . 

There is one other cusp to consider in (|5.11|) and figure ^ where p = 0. This 
cusp will somehow be smoothed on small scales affecting short distance physics. This 
will definitely alter the core energy of the string ( p. 231 ). It will not, however, affect 



the long distance behaviour of the string, which is typically the most important 
contribution and which is ultimately the source of the large A''^ dependence. 

5.5 Large N^. limit: summary 

As we have suggested earlier, unless there are peculiar degrees of freedom that greatly 
increase the domain wall tension, all of our results come under theoretical control 
in the large Nc limit Nc oo. The motivation for this comes from figure |^. In 
the large limit, the central barrier becomes extremely high as it is of order A^^ 
while the trough becomes flat. Thus, we approach the picture of figure |I] where the 
probability of nucleation becomes zero. In this limit, the domain walls become stable. 
To demonstrate this, we must show that, in this limit, the probability of creating a 
hole in the wall falls to zero by verifying the assertions (|5.4| ) and (|5.6|) . 
First, consider the dependence of the string tension a ( |5.23|) : 



a ~ "core + -^il log "^7,'^ ■ (5-33) 



The last term provides an A'c dependence of at least A'^logA'c. The core term can 
possible increase the tension, however, as we wish to demonstrated that the decay 
rate is suppressed, we must assume the worst possible case and take the lowest 
possible bound for a. Thus, we consider only the last term and set 

Ar,< Ar,logAr,<a. (5.34) 

This bound is quite robust. 

Now consider the three contributions to the domain wall tensions ( 5.27]) and 



(T^ ~ jlm^ ~ Nc^frfTq, (5.35) 
-}ElNl^^,, (5.36) 



'I 

~ m^fl/Nl ~ 1 . (5.37) 

To preserve the qualitative effects of the chiral limit, we first take this limit. Thus, 
despite of the fact that the vr meson cloud is of a much larger size than the r]' 
transition, its contribution to a is much smaller m^r ~ ^ in both the large, 

and the physical A'^c limits. Furthermore, it is a widely accepted assumption that, 
in QCD, all physical masses are of order 1 in the large Nc limit. The rj' mass is of 
course an exception ( p.l8| ), but only in the chiral limit as we have taken here: at fixed 
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mq, even m^/ ~ 1 in the large limit. The heavy degree of freedom (, however, is 
related to gluonic physics, and should thus not be sensitive to the chiral limit. If this 
is the case, then ~ 1 and the t]' contribution dominates in the large Nc limit. 

The reason that we have drawn so much attention to the point that in QCD, 
all masses seem to be of order 1 is that similar domain walls were discussed in the 
context of supersymmetric QCD [Q. There, a similar analysis to that performed in 
Section resulted in a domain wall of tension a ~ 1 <^ analogous to our esti- 
mate ( ^.36 ). However, in supersymmetric theories, this is in direct contrast with the 



analytic BPS lower bound of order Nc. In it was subsequently conjectured that 



heavy degrees of freedom at the cusp give a contribution similar to (|5.37|) but that, 
in order to reconcile this result with the BPS limit, the relevant mass of the heavy 
degrees of freedom has a peculiar dependence on the number of colours: ~ A^c- 

In QCD, there is no analogue of the BPS bound which can guide us, and we 
are not convinced that any such degrees of freedom exist with ~ Nc. However, 
imagining that they might exist, we take as a worst case^^ 

\^c<^<Nc. (5.38) 



Combining the results of Section we have the following Nc dependences 

NclogNc, a ~ . (5.39a) 

Rc ~ 0V^log Nc , /i-' ~ x/iVe , (5.39b) 

5*0 ~ A^^ P ~ e"^' . (5.39c) 

Notice also that, although both the critical radius for nucleation Rc and the wall 
thickness increase, the ratio fiRc ~ log Nc^ 1 and the semiclassical approxima- 
tion used to derive the decay rate (|5.1|) becomes justified in the large Nc limit. 

Allowing for the possibility that the contribution cr^ of the heavy field has an Nc 
dependence as conjectured in supersymmetric QCD, the tunnelling rate should still 
be exponentially suppressed: 

a>Nc\ogNc, a^Nc, (5.40a) 

Rc>\ogNc, ^i-^r^^c, (5.40b) 

So >Nc, P< e-^^ (5.40c) 



-'^^Notice that, if we do not first take the chiral hmit, then the pion contribution to the domain 
wall provides a ^ Nc contribution that will dominate in the large Nc limit. Such behaviour is 
not related to some mass ttItt ~ Nc- Thus, the pion contribution gives a Nc behaviour when 
niq 7^ is held fixed. Might it be possible that the resolution to the problem of the BPS bound 
in the supersymmetric case lies in an analogous degree of freedom to the QCD pion field which is 
not directly connected to the U(l) anomaly? These degrees of freedom might have been neglected 
in the relevant effective theory, yet may supply the required Nc dependence without requiring a 
field with strange Nc mass dependence. While we cannot say more about this here, the possibility 
supports the widely believed assumption that all masses in QCD are of order 1 in the large Nc 
limit. 
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In this case, our use of the semi-classical approximation — with independent cal- 
culations of the string and domain wall tensions — is no longer justified, however, 
the exponential contribution will still remain. The only way to ruin the qualitative 
picture would be to argue that a minimal contribution to the wall tension a exists 
that is > nI^'^ which is required to support the rj' transition = 27r. At present, 
we see no evidence or justification for such a contribution. 

5.6 iVc = 3 

We have argued that, in the large A*";, limit while preserving approximate chiral 
symmetry, QCD domain walls are stable to leading order and that, to leading order, 
they can decay through a tunnelling mechanism. We now provide some estimates 
of the magnitude of these effects extrapolating back to the realistic limit of N^. = 3. 
In this limit, we assume that the most relevant contribution to the domain wall 
tension is from the 7]': the pion contribution is suppressed by the ratio rriT^/mr^i 
and numerically contributes only 10% of the tension. Neglecting the contribution 
from the heavy field C, is more difficult to justify: numerically, both this and the rj' 
contributions are likely important, but it is reasonable to neglect them for an order 
of magnitude estimate. Thus we take 

cr = — 1 — cos — — . (5-41) 



V2 V 2Ar, 

In this formula, the \/2 in denominator should be replaced by ^jNj for an arbitrary 
Nf] however in all numerical estimates, we shall use Nf = 2 which we believe is 
very good approximation in the limit m„ ~ rrid <^ rris as equations (|2.1C ) and ( p.ll| ) 



suggest. Besides that, for Nc > 3, one can approximate 1 — cos[7r/ {2Nc)] — vr^/ {8N^) 
such that Equation ( |5.41| ) takes the simple form 

a = ■'ZJ 5.42 

which will be used for our numerical estimates. 

As an estimate for the string tension, we make the following estimate based on 
dimensional grounds: 

ar^V2E. (5.43) 
This estimate has the same Nc dependence as (|5.23| ).^''' 



^^The magnitude for a ~ V2E ~ (0.28GeV)^ should not be considered as a strong overestimation. 
Indeed, if one considers the tt meson string ||55|] which should be much softer (and therefore, would 
possess much smaller ajr) one finds, nevertheless, that a^r — tt/^ is very close numerically to this 
estimate for the rj' string tension. 
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For numerical estimates, we set Nc = 3, Nf = 2, and use: 



a ~ V2E ~ (0.29 GeV)^ (5.44a) 

a ~ ^^y^ ~ W.26 GeVf ~ (0.21 GeV)^ (5.44b) 
2v 2A'c 

= - ~ ^, (5.44c) 

In 



Although the estimate ( |5.44d| ) must be treated as very rough, it is nevertheless 



quite remarkable: In spite of the fact that all parameters in our problem are of 
order Aqcd, the classical action 5*0 may still be numerically large, and thus the 
corresponding tunnelling probability ( p.l|) might be quite small. We do not see any 
simple explanation for this phenomenon except for the fact that expressions (|5] 



and ( 5.44d ) for 5*0 contains a huge numerical factor 167r/3 of purely geometrical 
origin. In addition, since 5*0 ~ N^, one expects an additional enhancement in 5*0, 
even for A^c = 3. 

At this point we must address the validity of extrapolating from the large 
approximation to the physical region of Nc = 3. Qualitatively we expect the physics 
to be the same — the potential ( ^.7|) displayed in figure ^ has the same qualitative 
form for Nc > 2: namely, the central peak is much higher than the troughs and clas- 
sically stable solutions are admitted. As long as AVi > 0, the domain wall solutions 
remain classically stable, susceptible only to nucleation by string: the qualitative 
picture thus remains the same. Including heavier degrees of freedom, while certainly 
affecting the numerical results, should not modify this qualitative picture. 

Perhaps we should have included the effects of other light degrees of freedom 
such as the other pions, kaons etc. In the Nc = 3 case, these degrees of freedom 
are certainly relevant in general, but are not related to the physics which governs 
the U(l)^ QCD domain walls. We emphasize: QCD domain walls appear due to 
the discrete symmetry 6 —>■ 6 -\- 27t and the way that the singlet rj' meason uniquely 
interacts with the parameter 6 through the term {6—i log Det U). Thus, to satisfy the 
given boundary conditions, one expects the existence of a classical configuration for 
the 7]' field. The other light degrees of freedom exist as fluctuations in this classical 
background. Numerically one should account for these light particles scattering off of 
the domain walls, but for the qualitative estimates presented here, these light mesons 
are inconsequential. For the same reason, we do not expect the other mesons (p, uj 
etc.) to change the qualitative physical behaviour. The key point is the Mexican-hat 
shape of the potential ( |5.7| ) and that the potential of the central peak — where p = 
and the phase of the rj' phase is not well-defined — is much higher than any point 
in the trough defining the low-energy theory (|2.2|). This qualitative picture holds 
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in both the large Nc and Nc = 3 hmits. Low-energy fluctuations will not affect the 
height of the peak which makes the rj' domain walls classically stable. 

The 7]' mass serves to emphasize the validity of the large Nc limit in another 
way: The approximation of the chiral limit ( p.l6| ) remains valid up until A^^^ ~ 10. 
Thus, the physical case of A^^c = 3 should be related to the chiral large Nc limit that 
we have explored in this paper. 

One other interesting note: As we mentioned earlier, in QCD there is only one 
dimensional parameter, Aqcd? and it is thus generally believed that the semi-classical 
approximation in QCD cannot be parametrically justified. Nevertheless, Rc is nu- 
merically quite large — much larger than the width of the domain wall ( [4.13| ) 
which is set by m~,^. Therefore, the semi-classical approximation is somewhat 
justified a posteriori. 



5.7 Lifetime 

Now, we are prepared to make our last step and estimate the probability of creating 
a hole with radius Rc ~ 8Nc/ (tt^/^): 



rri C 




n 3 



e-^" X Det ~ nR'^Ei 




-So 



(5.45) 



We have estimated Det ~ i?4 dimensionally. Using the numerical values for Rc and 
5*0 — 130 given in ( |5.44| ) we arrive to the following final result: 

^ ~ lO^e'^^^GeV ~ lO^^^s"^ . (5.46) 

The most amazing result of the estimate (|5.46|) is the astonishingly small probability 
for the decay P ~ 10~^° GeV which one might naively expect to be on the GeV 
level. This small number leads to a very large life time for the domain walls, and 
consequently, makes them relevant to cosmology at the QCD scale. Of course, our 
estimation of 5*0 is not robust, and even small variation of parameters may drastically 
change our estimate (|5.46 ), making it much larger or much smaller. The point we 
wish to emphasize here is that it is at least possible for QCD domain walls to live 
long enough to have non-trivial physical effects. 

First of all, let us estimate an average size / of a domain wall before it col- 
lapses with an average lifetime of ~ //c. This corresponds the situation when the 
probability of the decay is close to one. If L is the Hubble size scale, L ~ 30 km, 
then 

3 



T 7ri?2 



P 

T ^ 



1. 



(5.47) 
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Formula (|5.47| ) implies that the average size I of the domain wall before it collapses 
could be as large as the Hubble size / ~ 10"^''^ ■ L ~ 0.2L. It also implies that on 
a Hubble scale domain wall L, the average number of holes which will be formed is 
approximately (n) ~ {L/lY ~ 25. Finally, the lifetime of a Hubble size domain wall 
is expected to be on the scale of 

Ti 7^~2-10~S, (5.48) 



which is macroscopically large! 

It is clear that all these phenomena are due to the astonishingly small num- 
ber (|5.46|) which makes the link between QCD and cosmology feasible. This small 



number is due to the tunnelling process rather than some special fine-tuning ar- 
rangements. As we mentioned above, we have not made any adjustments to the phe- 
nomenological parameters used in the estimates. Rather, we have used the standard 
set of parameters introduced in Equation ( p.2|) . To conclude the discussion of ( p.46| ) 
we would like to remind the reader that similar miracles related to tunnelling pro- 
cesses happens in physics quite often. For example, the difference in lifetime for U^^^ 
and Po^^^ under a decay is on the order of 10^° in spite of the fact that the "internal" 
physics of these nuclei, and all internal scales are very similar. 

We should be careful to qualify this result. First, it assumes that the effective 
Lagrangian ( p. 21) is a good description of the low energy physics of QCD including 
the physics of the rj' field. This assumption is well justified in the large Nc limit. 



The theory is not quantitatively justified as A^c 3 as we discussed in Section ^ 
Nevertheless, the qualitative features of the potential responsible for producing the 
domain walls persist for all A^'c > 3 as shown in figure |] so we expect the qualitative 
features to remain. Second, it assumes that the semi-classical calculation leading 
to ( p.45[ ) is justified. This is certainly not justified a-priori, but is somewhat justified 



a-posteriori. Third, we have assumed T = for this analysis. Thermal fluctuations 
may considerably reduce the lifetime of these walls. While this is a concern for 
cosmological domain walls, it should not be a concern at RHIC where the quenched 
approximation should be valid (see also footnote H). Finally, the interaction of the 
domain wall with nucleons can drastically change the numerical properties of the 
domain wall due to the strong interaction of all relevant fields. In particular, the 
domain walls may become much more stable in the presence of nucleons, and even 



account for the strong self-interacting dark matter discussed in p^ , |63| . We hope to 
return to these considerations in a separate publication. 

6. Conclusion and future directions 



The main results of this paper is expressed by the formulae (|4.12|) and ( ^.39|) where 



a new type of a quasi-stable QCD matter — the QCD domain wall (which could 
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have a macroscopically large size!) — is described. The quahtative picture for the 
QCD domain walls is exact in the chiral large Nc limit where the QCD domain walls 
are stable. To leading order, the walls decay as discussed, though the lifetime ( p.46| ) 
of the QCD domain walls is potentially very large. If this is the case, such objects 
might play an important role on the evolution of the early universe soon after the 
QCD phase transition. As a by-product, we also found a new type of axion domain 
wall (a^') with QCD structure. This may also have further cosmological applications 
if the axions exist. See for a first attempt in these directions. 

Unfortunately, we do not yet have a good microscopic picture of QCD in the low 
temperature, low density regime with small Nc = 3, thus we cannot further justify 
the effective theory ( p.2|) . This analysis, however, show that there is at least the hope 
of forming macroscopically large domain walls with QCD scale. Additional support 
for these objects comes from the high density regime we have better control of the 
physics. In this case, QCD domain walls nearly certainly exist Furthermore, 
if one accepts a conjecture on quark-hadron continuity at low temperature with re- 
spect to variations in the chemical potential fi then one can make the following 
argument: If for large /x domain walls exist, but for low /i they do not, then there 
should be some sort of phase transition as one lowers /i. Thus, the continuity conjec- 
ture supports the existence of quasi-stable QCD domain walls at lower densities, at 
least down to the densities of hypernuclear matter. Coupled with the fact that gluon 
and quark condensates do not vary much as one moves to the low density limit, we 
suspects that the qualitative picture holds even for zero density. 

It would be very exciting if QCD domain walls can be studied at RHIC (Rel- 
ativistic Heavy Ion Collider) [Q. Such a research would be the first experimental 
attempt to directly study the fundamental properties of the QCD vacuum structure. 
In this case the temperature would fall rapidly due to the rapid expansion and we 
expect that we can neglect thermal fluctuations as we have here. 

In any case, QCD domain walls are a very interesting feature of QCD and may 
live long enough to affect many areas of particle physics and astrophysics. We look 
forward to further study of these creatures! 
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